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We study the recently proposed effective-field theory for the phonon of an arbitrary nonrelativistic 
superfiuid. After computing the one-loop phonon self-energy, we obtain the low-temperature T 
contributions to the phonon dispersion law at low momentum and see that the real part of those 
can be parametrized as a thermal correction to the phonon velocity. Because the phonons are the 
quanta of the sound waves, at low momentum their velocity should agree with the speed of sound. 
We find that our results match at order T** In T with those predicted by Andreev and Khalatnikov for 
the speed of sound, derived from the superfiuid hydrodynamical equations and the phonon kinetic 
theory. We get also higher-order corrections of order T*, which are not reproduced pushing naively 
the kinetic theory computation. Finally, as an application, we consider the cold Fermi gas in the 
unitarity limit and find a universal expression for the low-T relative correction to the speed of sound 
for these systems. 

PACS numbers: 03.75.Kk;47.37.-|-q;11.10.-z 



I. INTRODUCTION 

Superfluidity is a phenomenon that occurs at low temperatures after the appearance of a quantum condensate, in 
either bosonic or fermionic systems [H, The condensate spontaneously breaks the global U{\) symmetry associated 
with particle number conservation of the system. In such a case Goldstone's theorem predicts the existence of low- 
energy modes that at sufficiently low momentum have a linear dispersion law and are essential to explain the property 
of superfluidity. We refer generically to these modes as superfiuid phonons, or phonons for simplicity. 

The phonons dominate the physics in the superfiuid at long wavelengths. At very low temperatures they also 
dominate the thermal corrections to the thermodynamical and hydrodynamical properties of the superfiuid. It is then 
important to have a precise knowledge of its self-interactions to know with accuracy the dynamics of these systems. 
Landau gave a successful phenomenological description of superfluidity and derived the main phonon self-interactions 
by imposing an ad hoc phonon Hamiltonian at leading order valid for "^He. Later on, it was shown that the energy 
spectrum for a weakly interacting Bose gas, which experiences a Bosc-Einstein condensation, could be derived |2|, |3|. 
The transport equation obeyed by the phonons can also be derived for the weakly interacting Bose gas [1] . 

Effective-field theory (EFT) techniques have been proposed and used with success in particle physics, although they 
have not been explored much in other branches of physics The techniques are especially suited for systems where 

there is a hierarchy of widely separated energy or momentum scales. Then, one performs an expansion in powers of 
energy or momentum, rather than expanding in a coupling constant. Power counting and symmetry considerations 
fix the form of the EFT to the accuracy one desires, regardless of whether the underlying system is weakly or strongly 
coupled. If the theory describes only the low-energy degrees of freedom, the dynamics of the short length physics 
is encoded in the value of the coupling constants of the effective theory. One remarkable successful example of the 
application of EFT techniques is chiral perturbation theory to describe the long-wavelength physics of the nuclear 
strong interactions 

EFT techniques have been recently applied to study superfiuid systems. The phonon Lagrangian at leading order 
in the momentum expansion and also the phonon Lagrangian at next-to-leading order have been explicitly derived, 
with an emphasis on applications to the superfiuid regime of the cold unitary Fermi gas [sl-flll [l3| . These are systems 
where the strength of the fermion interactions, measured in terms of the s-wave scattering length, is asymptotically 
large (sec Rcf. [l4| for a review). In this limit those systems exhibit conformal invariance, and they arc believed to 
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have universal properties [15|, meaning that their features do not depend on the detailed form of the interparticle 
potential. EFT techniques might be very conveniently applied to the superfluid phase of the cold unitary Fermi gas, 
as these are strongly coupled. 

This article is focused on the study of the recently proposed EFT to describe the phonons of a nonrelativistic 
superfluid and to extract from it their dispersion law. We discuss some subtleties that appear in the computations 
like, due to the fact that the phonons are massless degrees of freedom, the traditional EFT power counting gets mixed 
up when going close to the phonon on-shell limit. This problem shows up when studying the phonon dispersion law, 
as we discuss at length. We then focus on getting thermal corrections to the phonon dispersion law. 

We only consider the very low-temperature T regime of the superfluid, where the phonons are in a collisionless 
or nonhydrodynamical domain. In this situation, the phonons can be viewed as forming a bosonic gas, rather than 
a fluid. The hydrodynamical behavior of the superfluid in this situation simplifies, as opposed to what happens at 
higher T where the description of the system requires the two-fluid model of Landau. At T = 0, or very low T , there 
is only one fluid in the system, and one can talk about one single sound hydrodynamical velocity. 

At zero temperature, and very low momentum, the phonons obey 

lim Ep = CsP , (1) 
p— >o 

where Cg is the velocity of sound. This relation was suggested by Landau for the superfluid regime of ^He, as he 
realized that the elementary excitations in the system were the quanta of the sound waves Q . The relation has been 
proven for a Bose superfluid by Gavoret and Norieres in Ref. |16| for T ~ 0. Wc arc not aware of the existence of a 
general proof that the same relation should hold at vcry-low-temperature T. However, because the phonons are the 
quanta of the sound waves, it is expected that for very low momenta their properties should be the same as those of 
the hydrodynamical sound waves of the system. 

The speed of sound receives thermal corrections, first computed by Andreev and Khalatnikov [T^l in the collisionless 
regime of ^He, using both Landau's hydrodynamical equations and the phonon kinetic theory of Khalatnikov. This 
computation was believed to be valid for other superfluids as well. Thermal corrections to the phonon velocity 
have only been computed for the weakly interacting Bose gas [l^, [l^ , reaching agreement with the speed of sound 
correction in the terms of order T^lnT (see also Ref. [lO])- In this article we review the Andreev and Khalatnikov 
computation of the thermal corrections to the speed of sound and extend it to get the subleading corrections in 
T. We compare such a computation with the thermal corrections to the phonon velocity in the low-momentum limit, 
as derived from the EFT. We find agreement in the two quantities for the leading corrections that go as T'^lnT, but 
not for the subleading corrections of order T^. We discuss the possible origin of these discrepancies. 

This article is organized as follows. In Sec.|lT]we review the phonon EFT at leading order and next-to-leading order 
in a derivative expansion [l^. We present the analysis of the one- loop phonon self-energy with the leading-order 
Lagrangian in Sec. IIIII both at T = fSec. IIII X|) and finite T (Sec. IIII B)) . and discuss the need to go to the next-to- 
leading order for the computation of the phonon dispersion law in Sec. IIVI The phonon dispersion law derived from 
all the one-loop corrections is presented in Sec. |Vl In Sec. EH we review the computation of the thermal corrections 
to the speed of sound by Andreev and Khalatnikov and push it to get subleading corrections in T. In Sec. IVIII we 
concentrate on getting the thermal corrections to the speed of the phonon of the cold Fermi gas in the unitarity 
limit. Our conclusions are presented in Sec. I Villi Appendix [X] is devoted to showing some useful relations among 
the self-coupling constants of the phonon EFT, and in Appendix [B] wc show explicit details of the computation of the 
thermal corrections to the phonon self-energy. Throughout the article we use natural units, ?i = fc^ = 1. 



II. THE SUPERFLUID PHONON EFFECTIVE-FIELD THEORY AT LEADING AND 

NEXT-TO-LEADING ORDERS 

The superfluid phonon is the Goldstone mode associated with the spontaneous symmetry breaking of a U{1) 
symmetry, which corresponds to particle number conservation. EFT techniques can be used to write down the 
effective Lagrangian associated with the superfluid phonon. The effective Lagrangian is then presented as an expansion 
in derivatives of the Goldstone field, the terms of this expansion being restricted by symmetry considerations. The 
coefficients of the Lagrangian can be, in principle, computed from the microscopic theory through a standard matching 
procedure, and thus they depend on the short-range physics of the system under consideration. 

It has been known for a while that the leadin g-or der (LO) term Lagrangian of the Goldstone mode of a superfiuid 
system is entirely fixed by the equation of state |2ll423| . In some more recent publications |8|. |lO| it has been realized 
that at the lowest order in a derivative expansion the Lagrangian reads (see Sec. V of Ref. |1C| ) 



(2) 
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where P{fJ.o) and fiQ are the pressure and chemical potential, respectively, of the superfluid at T = 0, and 
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where tp is the phonon field, and m is the mass of the particles that condense. The form of X is dictated by imposing 
Galilean invariance to the EFT. In Ref. it was noted that after a Legendre transform, one can associate this 
formulation with Popov's formulation pTj . 

The final reason why the LO Lagrangian takes this particular form is that the effective action associated with the 
theory at its minimum for constant classical field configurations has to be equal to the pressure Q . This formulation 
turns out to be very advantageous, as it allows one to derive all the phonon properties at the lowest order in momentum 
based on the knowledge of the zero temperature pressure of the superfluid. In particular, one can easily get the 
phonon dispersion law and the form of the leading phonon self-interactions, and in principle their leading contribution 
to different physical processes. 

In order to see this, we expand the Lagrangian around P{p-o) in a Taylor series, and after the field redefinition 







(4) 



to have a canonically normalized kinetic term, one can write 



(5) 



We have neglected an irrelevant constant and a total time derivative term, which is only needed to study vortex 
configurations. 

All the coefficients that appear in Eq. ([5]) can be expressed in terms of different ratios of derivatives of the pressure. 
In particular, after using the thermodynamical relation dP = ^dfiQ , where po is the mass density at T = 0, one can 
check that at this order the phonon velocity is 
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that is, it can be identified with the speed of sound at T ~ 0, as is expected in the low-momentum limit. The 
dispersion law obtained from this Lagrangian, neglecting both quantum and thermal corrections, is exactly Ep = c^p. 
The remaining coefficients of Eq. ([5]) give account of the phonon self-couplings and are expressed as 
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In the following section we study the one-loop phonon self-energy, and we see that different combinations of the 
aforementioned coupling constants appear. For purposes of comparison with the kinetic theory approach, it is conve- 
nient to express these combinations in terms of derivatives of the speed of sound and of pq, rather than in terms of 
derivatives of the pressure. We present in Appendix [X] some of the relations that are used later on. 

The next-to-leading order (NLO) Lagrangian in the derivative expansion has been constructed by Son and Wingate 
in Ref. (see also Ref. [ll|) by coupling the system to both external gauge Aq and gravitational fields and demanding 
both coordinate and gauge invariance. For 9 = /iot — ip, it reads 



/:nlo = ^^X^^Xfl{X) + {V^ef f2{X) + mV^AofsiX) 



(8) 



where /i,/2 and /a are arbitrary functions and for simplicity we have assumed Aq = 0, as here we are not 
interested in the effect of an external potential in the phonon dynamics. For superfluid systems in the unitarity limit 
conformal invariance further restricts the form of the extra functions fi{X) [lo| and imposes /i(a*o) = ci{m/ po)^^'^ and 
/2(a*o) = C2(Ato/™)^^^i where Ci and C2 are dimensionless constants, which might be obtained from the microscopic 
theory by evaluating the static density and transverse response functions For the cold Fermi gas, these have 

been evaluated using the e expansion in Ref. |24l |. 
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FIG. 1: One-loop Feynman diagrams for the phonon self-energy, dubbed as bubble, tadpole and loUypop diagrams, respectively. 



When one works at the NLO, there are derivative corrections to the LO vertices, and also there are new vertices 
which are not present at the LO. But most important for the purposes of this article is that the phonon dispersion 
relation is now of the form 

Ep = c,p(l - 7/) + . . . (9) 
The parameter 7 is a function of the coefficients that appear in the Lagrangian at the NLO: 



+ ■ (10) 



The sign of the parameter 7 is very important in the phonon physics. If 7 < 0, then the process of one phonon decaying 
into two is kinematically allowed, but it is not allowed in the opposite case. This fact has several implications in the 
evaluation of the phonon damping, as wc show later on. 

The parameter 7 for the cold Fermi gas has been obtained from a phenomenological fit to numerical quantum Monte 
Carlo simulations in Ref. [2^ ; the authors of this reference fitted to the Monte Carlo data of Ref . [2^ . 

In principle, it would seem that the EFT at the LO described in Eq. 1^ should be enough to compute to leading- 
order quantum and thermal corrections. However, this is not always the case. When computing physical quantities 
close to the phonon on-shell limit, and due to the fact that phonons remain always massless degrees of freedom, one 
may detect sensitivity to the NLO dispersion law, even when computing the leading-order corrections. 

In the following sections we study the main corrections to the phonon dispersion law as derived from the EFT here 
described. We then see an explicit example where the problem discussed previously appears. 

III. ONE-LOOP PHONON SELF-ENERGY WITH THE LO LAGRANGIAN 

In this section wc compute the phonon self-energy at one-loop, as derived from the LO EFT described in the 
previous section. We study separately the nonthermal and the thermal corrections, which we denote as IIt=o{P) and 
IIt{P), respectively. As we will see, in the low-momentum limit, the first are very suppressed. 

The computation is done using the imaginary time formalism (IFT) (27j . That is, one first performs a rotation 
to Euclidean time. We denote with capital letters the four Euclidean momentum. Feynman rules for the phonon 
propagator and for vertices are straightforwardly derived from the effective Lagrangian equation [Eq.([5])]. We use 
dimensional regularization to deal with the ultraviolet (uv) divergences of the T = part of the diagrams, so that 
d = 4-e. 

The phonon propagator in the ITF with momentum = (ia;,i, q) = {qo, q) reads 

SiQ) ^ , (11) 

where = 27rnT, with n S Z, is a bosonic Matsubara frequency. At the LO, one has Eq = Csq, where q = |q|. 
There are three different diagrams that contribute to the one-loop self-energy. For external Euclidean momentum 

= {iuj,p) = {po,p), these are the bubble, the tadpole, and lollypop diagrams (see Fig. [T|). 
The bubble and the tadpole diagrams are expressed as 

nW(P) = -l%g^T j j^-^ [f,{P,Q)S{Q)S{P-Q)) , (12) 
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and 

n(*)(P) = -\T / ^^F,{P,Q)S{Q) , (13) 

respectively, where the vertex functions are 

Fi{p,Q) = [(po - 9o)pogo - % ((po - go)p • q + po(p • q- g^) + 9o(p^ - p • q))]^ , (14) 

F2{P,Q) = 12^2^2 _ 2r^;,.i(4pogoP • q + Po?" + 9oV) + 4r7A.2(2(p • q)' +pV) • (15) 
The lollypop diagram is expressed as 

n«(P) = 185^T /Ti^l^C^.O)^^^^'^^))^ E / A^(,^(/^,Q)5(/0), (16) 

where the vertex function is 

H{P, Q) = qopl - llg {qap^ + 2qoP ■ q) , (17) 
and the four-dimensional delta function is defined as 

S^^\Q) ^ ^S^,j2nfS^'\ci) , (18) 

where 6uj is a delta of Kronckcr. 

Tig 

A. One-loop self energy at T = 0. 

By dimensional analysis, it is easy to check that the one-loop corrections to the self-energy are proportional to 
p6 p at T = 0, and thus, they are very suppressed, particularly with respect to the thermal corrections we show 
later on. However, we still discuss how these small corrections might affect the phonon dispersion law at T = 0. 

In dimensional regularization the tadpole and lollypop diagrams are strictly zero. Thus, at T = we only need to 
consider the bubble diagram: 

J (271)^^(92 + c2g2)((cj-go)2+c2(p-q)2) 

Is it convenient to make use of Feynman parameters. Thus, we perform the change of variables R = Q — Px and 
rescale this last variable, doing rg = CsZq and r = z, and transform Eq. into 

tW /• d'^Z F,{P,Z) 



nVin(P) = — ^ / 7^ / dx ■ (20) 



We perform first the integral in the Z variable, using the well-known formula [26 

d'^z 1 _ 1 r(2-f) /'O^- ■ 



(277)'* (Z2 A)2 (47r)'*/2 r(2) \A 



(21) 



and similar ones with different powers of in the numerator. Here, as we are only interested in getting the lowest 
corrections to the phonon dispersion law, we approximate the numerator of Eq. ()20p . taking only the value for Z = 0, 
thus 

F,iP,0) = -.t2(1 - .t)2c.2(^2 ^ 3^^p2)2 (22) 



The point is that if one considers finite powers of Z in the numerator of Eq. (j20p and performs the Z integral, one 
gets more powers of A ^ (a;2 -|- c2p2) than those that appear in Eq. ([2T|) . Wc content ourselves with obtaining the 
lowest corrections to the dispersion law. In this limit, we also approximate Fi{P, 0) in a (aj2 ^ cj.p'^) expansion. Using 
the MS substraction scheme, and after analytical continuation uj = —ipo + 0+, we get 



rrib) I N '^sP ( pa dcs 



Po - c2p2 \ 47 



+ OW - <v') , (23) 



where we used the relations Eqs. (|A4[) and (jA5[) of Appendix El 



6 



B. Thermal corrections to the self-energy 



For the computation of the thermal component of the self-energy one can take d = 4, as it is uv finite. As we are 
only interested in obtaining the leading corrections to the dispersion law at low momentum, we only compute the 
self-energy in the limit pQ,p 0. In particular, we assume that these two scales are smaller that the scale set by the 
temperature. 

In the low-external-momentum limit several approximations can be performed. Due to the presence of the Bose- 
Einstcin distribution function, it is clear that the leading contribution to the integral comes from the region where 
q ^ T. In such a case, it is legitimate to approximate 



Fi{P, Q) « {-po{ql - Vgl^) - 277g9oP • q)^ 
in the numerator of the bubble diagram, and also 

dnB{Eg) 



El — E2 ~ CsP cos 6 , fi — f2 ~ Cs p cos ( 



dE„ 



(24) 
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We present in Appendix [B] the explicit details of the computation of the three Feynman diagrams, and here we only 
display the final result, valid in the po,p —i- limit: 
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(26) 



As clearly seen from Eq. (j23|) . for low external momentum, assuming that p <^ T, the T = self-energy is very 
suppressed with respect to the thermal part of the self-energy. 



IV. LEADING CORRECTIONS TO THE SELF-ENERGY COMING FROM THE NLO LAGRANGIAN 

In order to study both the thermal and the higher momentum corrections to the linear phonon dispersion law, 
knowledge of the phonon self-energy is necessary. Naively, one could think that computing the one-loop self-energy 
with the LO Lagrangian should be enough to get the leading corrections to Eq. (jlj. However, given the explicit form 
of both IIt{P) [see Eq. ([25)1 ] and IIt=o{P) [see Eq. we see that when studied close to the on-shell limit, these 

might be logarithmically sensitive to the NLO dispersion law as well. If the corrections provided by the one-loop LO 
physics to the phonon dispersion law are smaller than those associated with the value of the parameter 7. see Eq. ([9|). 
then one cannot neglect the NLO Lagrangian in the first corrections to the dispersion law. We will assume that this 
is the case, as the thermal corrections are governed by the ratio (T^/pocl), and the T = corrections are governed 
by the ratio {p'^/po), and these two ratios are very small in the low-temperature regime of the superfluid. In this 
situation, we must start our computation with the Lagrangian described by >Clo + ^NhO- 

In this section wc compute the leading thermal and nonthermal corrections of the theory described by £lo + 'Cnlo ■ 
Several new one-loop diagrams appear, but most of them will produce higher p or T corrections than those computed 
in the previous section. Those are the diagrams that arise from higher-momentum corrections to the different vertices 
of the LO physics, or new vertices with higher-momentum dependence. We do not compute these higher-order 
corrections here, but simply focus on getting the corrections that contribute at the same order as those computed 
with £lo- 

It is very easy to spot the one-loop contributions that give the leading one-loop corrections. Consider the bubble 
diagram, with vertices arising only from £lOi but with the propagator computed from £lo + -Cnlo- The phonon 
propagator takes the same form as in Eq. but now the energy is given by Eq. 

Let us first consider the thermal corrections to the self-energy. If we only consider the low-external-momentum 
limit, the analysis turns out to be very simular to the one we carried out in the preceding section. Carrying out the 
same steps and approximations, one only has to replace Eq. (|25p by 



El - E2 ~ CsP cose {1 ^ 3jq^) 



(27) 
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After an analysis very similar to the one we carried out before (see Appendix |B] for details) , we reach the following 
expression when the phonon is nearly on-shell, po ~ Csp{l — 7P^), 



nT(Po,P) 



15c5 pq 
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1 Po dcs 

2 Cs dpo ' c?^p'^ 
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A similar thing happens at T = 0. The bubble diagram is modified when considering the corrected phonon 



propagator. In particular, it now reads 
(b)JVLO/p^ 18g^ 
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where FsiP, Z,x) ^ -2[(z + a;p)**(l - a;) + (z - (1 - x)^fx\. This expression can be simplified by the following 
argument: the term proportional to 7 in the denominator will only be of the same order of magnitude as the rest of 
the denominator for very small Z^, but if this happens P ^ Z and this means that the approximation i^3(P, Z,x) ~ 
F3{P,0,x) can be made. If one considers also that the external particle is nearly on-shell,a;^ + c^p^ sa 2c^7p^, then 
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In the MS substraction scheme, and after analytical continuation to Minkowski space, we get 
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(30) 
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V. THE PHONON DISPERSION RELATION 



The phonon dispersion law is obtained by studying the poles of the resumed propagator, that is, after solving 

p2-p2^2_n(po,p)=0, (32) 

where n is the one-loop self-energy. Because the T = contributions to this quantity are very suppressed, we first 
only consider the temperature-dependent part. 

We divide Eq. ([5^ by p^c^ and realize that the real part of the equation only leads to thermal corrections to the 
phonon velocity. After writing po ~ (c, + 5v{T))p — ia{p, T), we reach 
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The leading thermal corrections we have obtained, that is, those that go as ^ T'' InT are the same as those obtained 
from the kinetic theory approach of Andrccv and Khalatnikov |T3l, which wc review in Sec. IVII 

For completeness, we also discuss the modifications to the dispersion law at T = that arise from the one-loop 
physics. If we assume po « CsP (l — 7?^^) + ^Ep — ia{p), then we obtain 

CsP^ 
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and 



a{p) = e(-7) 
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As for the damping terms, we note that if we work with the LO phonon Lagrangian, we obtain the same damping 
term, but without the step function of —7. For the case of a weakly interacting Bose gas, where the following 
property is fulfilled f^f^ = 5 (see, for example Ref. (sOj), we then get the results for the damping term first obtained 
m Refs. milllsl. However, when working to the NLO, with a dispersion law of the form expressed in Eq. the 
process described by a, which corresponds to one phonon decaying into two, is only kinematically allowed for negative 
values of 7. 

The thermal corrections to order In T of the phonon velocity of a weakly interacting Bosc gas have been computed 
in Ref. |20| Jsee Eq. (54)]. We note that we agree with the value of the coefficient in front of the logarithm. However, 
Giorgini [20[ computed only at the LO and thus does not show any dependence of the result on 7, while our computation 
proceeds at the NLO and shows the sensitivity of these logarithmic corrections to the nonlinear corrections to the 
phonon dispersion law, which cannot be neglected. 



VI. THERMAL CORRECTIONS TO THE SPEED OF SOUND IN A COLD SUPERFLUID FROM 

KINETIC THEORY 

The thermal corrections to the speed of sound in a ^He superfluid at low T were first computed by Andreev and 
Khalatnikov [13] in a collisionless regime. For completeness, wc review such a computation here (sec also Ref. 0]) 
and extend it to get also the subleading thermal corrections. 

We consider Landau's two fluid equations for the superfluid [ll, Q to study the evaluation of the dynamics of 
sound waves. For our purposes, it is enough to concentrate only on the continuity equation and the equation for the 
superfluid velocity Vgi: 

dtp + divj =. , (37) 



dtv, + V (^Ai + Y j = , (38) 

where p is the mass density, j is the current, and p is the chemical potential. 

At zero temperature only the pure superfluid component exists. At finite temperature, a normal fluid component has 
to be considered as well, which takes into account the contribution to the hydrodynamics of the different quasiparticles. 
However, at very low temperatures only the massless quasiparticles of the superfluid, the phonons, are thermally 
excited and contribute to the thermodynamical and hydrodynamical properties of the system. 

In the very-low-T regime the mean free path of the superfluid phonons is very large [4] , and one can consider that 
they are in a collisionless regime or, in other words, that they are not in a hydrodynamical regime. In this case the 
thermal phonons can be viewed as forming a bosonic gas, and thus they can be described with a Boltzmann equation. 

In the collisionless regime the phonon distribution function n obeys the kinetic equation Q 

dn ^ dn dH dn dH ^ 
dt dr dp dp dr ' 

where H ^ Ep + p ■ Vs is the phonon Hamiltonian. The energy of the phonon can be taken as in Eq. © if we move 
from the strict low-momentum limit. The phonon transport equation for the weakly interacting Bose gas has been 
derived in Ref. g (see also Ref. [|). 

The aforementioned transport equation has been explicitly derived for a weakly interacting Bose gas, starting from 
the microscopic description of the gas Q . We are not aware of a similar derivation in the literature for other superfluid 
systems. 

The phonons give contributions to the current and chemical potential of the system, which are expressed as 

f f 9Eg 

j^^"^+y (2^^"' '^^'^°+y (2^^"' 

where /ig is the chemical potential at T = 0. 

One then studies deviations from the equilibrium solutions, proportional to g-^'^i+^k r^ that 

n ^ n + Sn , p — p + Sp, (41) 

where the quantities with the bar refer to the equilibrium values. For the unperturbed value of the phonon distribution 
function we take the Bose-Einstein distribution function, n = nsiEp). 

We also assume that we work in the frame where Vg = 0, but there is also a fluctuation in the superfluid velocity, 5vs. 
In Ref. Q , it is mentioned that one should also consider the functional dependence of the phonon energy on the density 
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of the excitations, with a term of the form 5Ep = J , f{p,p')5n{p'), where / is a function to be determined. This 
term was neglected by Andreev and Khalatnikov jl7{" under the assumption that it gives very subleading corrections 
to the final result. For the time being, wc make the same assumption. However, in order to check the accuracy of 
this assumption, one might need to evaluate the form of the function /, starting from the microscopic description of 
the superfluid. 

After linearizing the kinetic equation (j39p . one finds the value of the fluctuation in terms of 6p and (Jv^: 



Sn ~ — 



dni 



Vph 



dEp w — k • Vph 



dEp 
dp 



(5p + p • (5v. 



(42) 



where the phonon velocity is given by Vph — -g^- Consequently, one can compute the fluctuations in both the current 
and the chemical potential due to these phonon fluctuations using Eqs. (|40|) . After linearizing also the hydrodynamical 
equations, Eqs. p7p and psp. one gets a homogeneous system of linear equations for 6p and (5vs. Keeping only the 
leading and subleading thermal corrections, and assuming Wph ~ Cg, one flnds 



UJ 



p^dp duB 
27r2 Ie^ 

dp 



pd'^Cs 



'4 9p2 





)■-} 


(5- 













- 2 



iVph 



ph 




dp fc^Wph 
dcos 9 



^0^ 



(43) 



and retarded boundary conditions have been assumed, with the +iO prescription. 

Equation ()43p might be interpreted as the dispersion law for the hydrodynamical wave. Its solution contains both 
real and imaginary parts, and the first is seen to correct only the value of the speed of sound. Thus, if we define 
a; (cg + Scs(T))k — ia{T, k), where 5cs is the correction to the sound velocity due to the thermal fiuctuations and 
a is the sound decay coefficient, or attenuation factor, one then gets 



5c,{T) 



Ic. 
2 p 



p dp dn B 
27r2 dK 



1 P^ d^Cs 

4 Cs dp^ 



P dcs 
Cs dp 



ph 



1 , UJ — kv, 
- log 

2 UJ + fcwph 



(44) 



Andreev and Khalatnikov [T3| solved the dispersion relation, Eq. (|44p . keeping only the leading terms in T. They 
noticed that the solution was sensitive to the form of the phonon dispersion law. Eq. and the parameter 7. The 
equation can also be solved keeping subleading terms in T. One finds 



Sc4T) 



TT 

30^ 



T 



1 



p dcs 
Cs dp 



lOE 



27|7|T2 




We also find a nonvanishing attenuation factor if the parameter 7 is negative: 



„(r,.,.e(-,)--j|i + ^^ 



P dc, 
Cs dp 



(45) 



(46) 



where O is the step function. Andreev and Khalatnikov [l^l did not consider such a term, as they assumed that 7 
was positive for '*He. 

We agree with Andreev and Khalatnikov in the correction of order T'^lnT to the speed of sound, except 
for a factor 2 in the argument of the logarithm, as already pointed out in Ref. [l^: which is probably due to the 
approximation that these authors used to solve the integral of Eq. (j44[) . Here instead the corresponding integral was 
computed exactly. 

Because at low T one can approximate p w po: we see that the corrections of order T'^lnT are the same as those 
obtained with the EFT for the phonon velocity [see Eq. p3p ]. As for the corrections of order T^, there is a clear 
discrepancy between the two results. However the discrepancy involves the same kind of thermodynamical derivatives, 
an indication that the two approaches must be very similar. We believe that this discrepancy might require one to 
have a rigorous derivation of the phonon kinetic theory approach for every superfluid system, so as to derive all 
possible sources of the subleading thermal corrections to different physical quantities. 
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VII. APPLICATIONS TO THE COLD UNITARY FERMI GAS 



111 the unitarity limit, the thermodynaniie properties of the cold Fermi gas can be determined up to some dimen- 
sionless constants (T5| . At zero temperature, and due to the absence of any internal scale, dimensional analysis fixes 
the form of the pressure as being proportional to that of a free system, 

P = com3/Vr, (47) 
where cq is a dimensionlcss and universal constant. This parameter can be expressed as 

25/2 

^0 = - (48) 

where ^ is the universal constant that fixes the relation between the chemical potential and the Fermi energy fiQ — S,Ep. 
Experiments with cold trapped fermionic atoms (33l | find £ ~ 0.32 — 0.44, a result that is in agreement with Quantum 
Monte Carlo calculations at vanishing temperature [13, Hg. 

At sufficiently low temperatures, the system is a superfluid phase. The temperature for the superfluid phase 
transition has been determined also with Monte Carlo numerical simulations, and then one finds Tc ^ Q.23Ef ~ ^^^Mo 

At zero temperature one can easily deduce the density from the value of the pressure, as well as the value of the 
speed of sound, which is = 2^o/3m. In this case it is easy to check 

Cs dpo 3 ' Cs dpi 9 

and that the combination 

only depends on fia and on £, but not on the value of the mass of the specific fcrmion species one is considering. 

The values of £ and 7 can be extracted from a fit to numerical Monte Carlo simulations [25| . Following the notation 
of Ref . [2^ , we can parametrize 

1= (51) 

c2 Sm^ci 'S2pl ' > 

and A is also an adimensional universal constant. The fit to numerical simulations that was performed in Ref. [25j 
gives the result £ = 0.455 and A = 0.13. 

We can express the relative thermal correction of the phonon velocity as 



^^-^^^f^V(43 + 161ogf^^) 
Cs 160 \ ^\ 32 ^il J 



(52) 



which depends only on the ratio and the constants £ and A. It is important to stress the universal character of 
this thermal correction for all cold Fermi gases at unitarity. 

In Fig. m wc plot Eq. ((5^ . using the numerical the values of £ and A obtained in Ref. [25|, up to temperatures 
below Tc, the critical temperature for the transition to the normal phase. As we see from the plot, the corrections 
of order compcnsantc those of order In T and make the thermal correction to the phonon velocity really tiny 
at still relatively low temperatures. At sufficiently low temperatures, one may expect that these corrections arc the 
same as those for the speed of sound, as the phonons are the only thermal excited states. At temperatures close to Tc 
one may expect a discrepancy between these two quantities, as other quasiparticles (for example, the fcrmion pairs) 
may contribute to the corrections to the speed of sound. 

The sound attenuation factor is also given by 

a(P,T) = — J . (53 

135pocf 

At larger temperatures, when phonons are in a hydrodynamical regime, the expression for the sound attenuation is 
different, as it is expressed in terms of the different transport coefficients of the system 0. For cold Fermi gases the 
sound attenuation in the phonon hydrodynamical regime has been recently computed in Ref.jl^. It is not clear at 
which values of T one could consider that the phonons are in a hydrodynamical regime or not, as this would depend 
on the ratio between the phonon mean free path and the macroscopic volume of the system. 
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FIG. 2: Relative thermal correction to the phonon velocity for a cold Fermi gas in the unitarity limit. With the dashed line 
we plot the corrections of order In T, which agree with that of the speed of sound predicted with kinetic theory. With the 
solid line, we plot the entire correction that includes the subleading terms of order T^, given in Eq. (|52|) . We use the values 
^ = 0.455 and A = 0.13 obtained in Ref. 

VIII. CONCLUSIONS 

In this article we have studied the phonon dispersion law in an arbitrary superfluid system as derived from the 
phonon EFT in a collisionless or, equivalently, very-low-T regime. We have concentrated our efforts on considering 
how the speed of the phonon is modified at low temperature. Because the phonons arc the quanta of the sound waves 
in the superfluid, when looking into this modification in the low-momentum limit, such a thermal correction Sv(T) 
should be the same as the one corresponding to the speed of sound, 5cs{T), as in this regime the phonons are the 
only excited states. 

We have compared the results of Sv(T) with the old computation of Scs{T) due to Andreev and Khalatnikov \vi\ . 
which is based on studying the superfluid hydrodynamical equations, together with a kinetic transport equation for 
the phonons. The last was deduced based on very general principles, assuming the ad hoc phonon Hamiltonian of 
Landau. It has only been derived from the microscopic theory for the superfluid phase of a weakly interacting Bose 
gas. For very low T, our results agree with those of Andreev and Khalatnikov p7j in the corrections of order T*^ InT. 
However, when we push the kinetic theory approach to higher-order corrections in T, we find a different result than 
that predicted with EFT techniques. We believe that the origin of this discrepancy might require a modification 
in the standard kinetic theory approach. Possible corrections might be due to the change in the phonon spectrum 
in the presence of fluctuations of the density of the phonon excitations, a term already discussed in Ref. 0. It 
is, however, not clear how to derive such corrections in the kinetic theory approach without studying the specific 
underlying microscopic system under consideration. On the other hand, the EFT techniques we used are based on 
simple symmetry and power counting arguments, which are universally valid for all superfluid systems. 

Apart from demonstrating the efficiency of the EFT techniques, our ultimate goal is considering a well-defined 
framework to study the phonon contributions to the thermodynamical and hydrodynamical properties of the su- 
perfluid, especially those which are not weakly coupled, such as the Fermi liquid close to the unitarity limit. Very 
recently, different transport coefficients for the superfluid regime of these systems have been computed IstHsq} . The 
computations were done using the kinetic theory described by Khalatnikov [2[ and used here. While all computations 
were done to the leading order in a temperature expansion, it might be interesting to see whether modifications 
to the kinetic approach are needed if one wants to study subleading corrections in T of the transport coefficients. 
Alternatively, one could compute the different phonon contributions to the transport coefficients using the phonon 
EFT. However, as known for other systems such as quark-gluon plasma [i^, the consideration of collision terms in 
a field theory treatment would require one to study the theory at the two-loop order. We leave all these important 
issues for future studies. 



Appendix A: Some useful relations 

In order to compare the results of the phonon dispersion law with the results obtained for the dynamics of the 
sound waves, it is convenient to express differently the phonon coupling constants of the LO Lagrangian that appear 
in Eqs. ([7]). While in Eqs. ^ those arc expressed in terms of ratios of derivatives of the pressure with respect to the 
chemical potential, it will turn out to be more convenient to express them as derivatives of the speed of sound with 
respect to the density. 

Like in the main part of this article, in this appendix we denote quantities evaluated at zero temperature with the 
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subscript 0. From the definition of tlie speed of sound c^, one can check that at zero temperature one can write 

iM^i 1 d'P _ (Al) 

2 9/io 2m 2771 / Q2p\^ dfil 

Also from the same definition and that of the mass density po, one has 

^ = m^=™^^. (A2) 

Using the previous two equations, one can immediately check that 

Po dcs I _ 1 |£ d^P 1 cl 
Cs dpo^ 2 2 (a2p\'^ diil ~ 2 2rig ' 



(A3) 



where rjg is the phonon self-coupling defined in Eqs. ([7]). From this expression we deduce the following relation 



2 



4r;2 V ^' \csdpo 

From the definition of the coupling constants, see Eqs. (O, it is easy to infer 



and 



5' = (A5) 



A.., = 1^ . (A6) 

Also from Eqs. ^ and (|A3[) one can check that 

£l!?M = 2i = 2-4^|^. (A7) 
Another useful relation that is needed in the one-loop computations performed in this article is 

'7A,2 3 \^ CsdpQ cj \dpaj dp^ J 

which can be deduced after some algebra, from the definitions Eqs. ([7]), and after expressing ^-f in terms of derivatives 
of the speed of sound and of po ■ 



Appendix B: Computation of the thermal one-loop self-energy 

We present in this appendix the explicit details of how to obtain the thermal component of the one-loop phonon 
self-energy in the limit of low external momentum. 

We first consider the bubble diagram. After performing the sum of Matsubara frequencies one finds [29l | 

ui^HP) ^ ISg^ y: f^F.iP,Q) ( ^l±fi^lEll±nfl^] , (Bl) 

where E\ = Eq and E2 = ^'|p-q|, and 
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We use the same notational conventions of Ref. (27[, so that the equihbrium Bose-Einstein distribution function is 
nsiE) = f{\E\). Further one can see that f{—E) = —(1 + nB{E)). We evaluate separately the pieces of Eq. (|Bip 
where si — ~S2 and where si = S2, taking into account only the thermal contributions to the integrals. For the pieces 
where si = —S2, we have 



(fiq dns 1 



(27r)3 dEq E^ iu) — CspcosO 



[Fi{P,Q)+F,{P,-Q)] 



<3o=-E, 



(B3) 



where we have used the approximations given in Eq. (|25p . and we also consider that in the low- momentum external 
limit one can use Eq. (|24p . Now we proceed in exactly the same way with the pieces of the bubble diagram, Eq. (|Bip , 
when si = S2, considering only the thermal part of the diagram. In the limit of low external momentum wc get 



ui%{p) 



9.g2 f d^q UBjEq) 
(27r)3 



[Fi(P,Q)+Fi(P,-Q)] 



qo=Eg 



(B4) 



At this point, we also perform the analytical continuation to Minkowski space with retarded boundary conditions, 

loj — > Po + ^O"*", and thus we get 



ni%(po,p) 



1 



d^Q duf 



cos 



'^pUvi - '^Vgci + ct) + Spocspiv' - Vgci) cose + ^ifgclp" ■ 



^t^IpoJ (2^)3 dE -^-cos0 + iO+ 



:os2 e] 
(B5) 

where we used Eq. (jA5[) . 

It is then realized that this part of the diagram has both real and imaginary parts, the last one being related to 
the damping of the phonon. Using the relation Eq. (|A4p of Appendix I3 it can be rewritten as 



nfkpo,p) 



Po 



d^q duB 2 
(27r)3 dEq"^^ 



Po 

CsP 



Pa dcs 
Cs dpo 



Po dcs 
Cs dpo 



Pa 



CsP 



(B6) 



As for Eq. (|B4[) . using the the relation Eq. (|A5[) . together with Eq. (jA3p . it can be expressed as 



nC') (r. r.^ ~ ^^^^ / _f^dnB 2 

^h,T\POiP) :— J j2nf~^^ 



12 4 \csp) c2 \dpo 



Po J [2tt)^ dE 

The frequency sum in the tadpole diagram Eq. (jl3|) can be done using standard techniques and gives 

J [^^r qo=sEg 2Eq 



(B7) 



(B8) 



After analytical continuation to Minkowski space, and using the relations Eqs. (jA6p . (|A7[) . and (jASp of Appendix [XI 
it transforms into 



4*^(po,p) 



cjp^ 



d'^q duB o 
-9 



J_ Po dcs 
12 4cs dpo 



Po 

CsP 



po J (27r)3 dEq ■ 
The lollypop diagram is computed in a similar way, and it gives 



2 / o_ \ 2 


Acs dpo ' 4c2 yOpo^ 



3po dcs bpl ( dc. 



pI d'^Cs 
4c, dpi 



n^T\po,P) 



n^T^p^ Po dcs 
30c3po Cs dpo 



Po 

CsP 



2jOo dcs 
Cs dpo 



(B9) 



(BIO) 



where we have performed first the q integral in Eq. (|16p and then the sum over Matsubara frequencies. Inverting 
the order of these two last operations gives an ill-defined result. We have verified our computation by analyzing 
the diagram in configuration space, using a mixed representation of the phonon propagators, where no ambiguity is 
present. 

Adding the different contributions obtained so far. 



iit{po, p) = TirApo, p) + n^"^(po, p) + ii^t\po, p) + (po, p) 



i['^(po,p) + ng 

and computing the momentum integrals, one reaches the final result given in Eq. (|26p 



(Bll) 
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